Development of a Droplet Digital PCR to Monitor SARS-CoV-2 Omicron Variant BA.2 in Wastewater Samples.

Wastewater-based surveillance can be used as a complementary method to other SARS-CoV-2 surveillance systems. It allows the emergence and spread of infections and SARS-CoV-2 variants to be monitored in time and place. This study presents an RT-ddPCR method that targets the T19I amino acid mutation in the spike protein of the SARS-CoV-2 genomes, which is specific to the BA.2 variant (omicron). The T19I assay was evaluated both in silico and in vitro for its inclusivity, sensitivity, and specificity. Moreover, wastewater samples were used as a proof of concept to monitor and quantify the emergence of the BA.2 variant from January until May 2022 in the Brussels-Capital Region which covers a population of more than 1.2 million inhabitants. The in silico analysis showed that more than 99% of the BA.2 genomes could be characterized using the T19I assay. Subsequently, the sensitivity and specificity of the T19I assay were successfully experimentally evaluated. Thanks to our specific method design, the positive signal from the mutant probe and wild-type probe of the T19I assay was measured and the proportion of genomes with the T19I mutation, characteristic of the BA.2 mutant, compared to the entire SARS-CoV-2 population was calculated. The applicability of the proposed RT-ddPCR method was evaluated to monitor and quantify the emergence of the BA.2 variant over time. To validate this assay as a proof of concept, the measurement of the proportion of a specific circulating variant with genomes containing the T19I mutation in comparison to the total viral population was carried out in wastewater samples from wastewater treatment plants in the Brussels-Capital Region in the winter and spring of 2022. This emergence and proportional increase in BA.2 genomes correspond to what was observed in the surveillance using respiratory samples; however, the emergence was observed slightly earlier, which suggests that wastewater sampling could be an early warning system and could be an interesting alternative to extensive human testing.

A general approach to identify low-frequency variants within influenza samples collected during routine surveillance.

Influenza viruses exhibit considerable diversity between hosts. Additionally, different quasispecies can be found within the same host. High-throughput sequencing technologies can be used to sequence a patient-derived virus population at sufficient depths to identify low-frequency variants (LFV) present in a quasispecies, but many challenges remain for reliable LFV detection because of experimental errors introduced during sample preparation and sequencing. High genomic copy numbers and extensive sequencing depths are required to differentiate false positive from real LFV, especially at low allelic frequencies (AFs). This study proposes a general approach for identifying LFV in patient-derived samples obtained during routine surveillance. Firstly, validated thresholds were determined for LFV detection, whilst balancing both the cost and feasibility of reliable LFV detection in clinical samples. Using a genetically well-defined population of influenza A viruses, thresholds of at least 104 genomes per microlitre and AF of ≥5 % were established as detection limits. Secondly, a subset of 59 retained influenza A (H3N2) samples from the 2016-2017 Belgian influenza season was composed. Thirdly, as a proof of concept for the added value of LFV for routine influenza monitoring, potential associations between patient data and whole genome sequencing data were investigated. A significant association was found between a high prevalence of LFV and disease severity. This study provides a general methodology for influenza LFV detection, which can also be adopted by other national influenza reference centres and for other viruses such as SARS-CoV-2. Additionally, this study suggests that the current relevance of LFV for routine influenza surveillance programmes might be undervalued.

Optimization and Application of a Multiplex Digital PCR Assay for the Detection of SARS-CoV-2 Variants of Concern in Belgian Influent Wastewater.

Since the beginning of the COVID-19 pandemic, the wastewater-based epidemiology (WBE) of SARS-CoV-2 has been used as a complementary indicator to follow up on the trends in the COVID-19 spread in Belgium and in many other countries. To further develop the use of WBE, a multiplex digital polymerase chain reaction (dPCR) assay was optimized, validated and applied for the measurement of emerging SARS-CoV-2 variants of concern (VOC) in influent wastewater (IWW) samples. Key mutations were targeted in the different VOC strains, including SΔ69/70 deletion, N501Y, SΔ241 and SΔ157. The presented bioanalytical method was able to distinguish between SARS-CoV-2 RNA originating from the wild-type and B.1.1.7, B.1.351 and B.1.617.2 variants. The dPCR assay proved to be sensitive enough to detect low concentrations of SARS-CoV-2 RNA in IWW since the limit of detection of the different targets ranged between 0.3 and 2.9 copies/µL. This developed WBE approach was applied to IWW samples originating from different Belgian locations and was able to monitor spatio-temporal changes in the presence of targeted VOC strains in the investigated communities. The present dPCR assay developments were realized to bring added-value to the current national WBE of COVID-19 by also having the spatio-temporal proportions of the VoC in presence in the wastewaters.

Strategy to Develop and Evaluate a Multiplex RT-ddPCR in Response to SARS-CoV-2 Genomic Evolution.

The worldwide emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) since 2019 has highlighted the importance of rapid and reliable diagnostic testing to prevent and control the viral transmission. However, inaccurate results may occur due to false negatives (FN) caused by polymorphisms or point mutations related to the virus evolution and compromise the accuracy of the diagnostic tests. Therefore, PCR-based SARS-CoV-2 diagnostics should be evaluated and evolve together with the rapidly increasing number of new variants appearing around the world. However, even by using a large collection of samples, laboratories are not able to test a representative collection of samples that deals with the same level of diversity that is continuously evolving worldwide. In the present study, we proposed a methodology based on an in silico and in vitro analysis. First, we used all information offered by available whole-genome sequencing data for SARS-CoV-2 for the selection of the two PCR assays targeting two different regions in the genome, and to monitor the possible impact of virus evolution on the specificity of the primers and probes of the PCR assays during and after the development of the assays. Besides this first essential in silico evaluation, a minimal set of testing was proposed to generate experimental evidence on the method performance, such as specificity, sensitivity and applicability. Therefore, a duplex reverse-transcription droplet digital PCR (RT-ddPCR) method was evaluated in silico by using 154 489 whole-genome sequences of SARS-CoV-2 strains that were representative for the circulating strains around the world. The RT-ddPCR platform was selected as it presented several advantages to detect and quantify SARS-CoV-2 RNA in clinical samples and wastewater. Next, the assays were successfully experimentally evaluated for their sensitivity and specificity. A preliminary evaluation of the applicability of the developed method was performed using both clinical and wastewater samples.

Strategy and Performance Evaluation of Low-Frequency Variant Calling for SARS-CoV-2 Using Targeted Deep Illumina Sequencing.

The ongoing COVID-19 pandemic, caused by SARS-CoV-2, constitutes a tremendous global health issue. Continuous monitoring of the virus has become a cornerstone to make rational decisions on implementing societal and sanitary measures to curtail the virus spread. Additionally, emerging SARS-CoV-2 variants have increased the need for genomic surveillance to detect particular strains because of their potentially increased transmissibility, pathogenicity and immune escape. Targeted SARS-CoV-2 sequencing of diagnostic and wastewater samples has been explored as an epidemiological surveillance method for the competent authorities. Currently, only the consensus genome sequence of the most abundant strain is taken into consideration for analysis, but multiple variant strains are now circulating in the population. Consequently, in diagnostic samples, potential co-infection(s) by several different variants can occur or quasispecies can develop during an infection in an individual. In wastewater samples, multiple variant strains will often be simultaneously present. Currently, quality criteria are mainly available for constructing the consensus genome sequence, and some guidelines exist for the detection of co-infections and quasispecies in diagnostic samples. The performance of detection and quantification of low-frequency variants using whole genome sequencing (WGS) of SARS-CoV-2 remains largely unknown. Here, we evaluated the detection and quantification of mutations present at low abundances using the mutations defining the SARS-CoV-2 lineage B.1.1.7 (alpha variant) as a case study. Real sequencing data were in silico modified by introducing mutations of interest into raw wild-type sequencing data, or by mixing wild-type and mutant raw sequencing data, to construct mixed samples subjected to WGS using a tiling amplicon-based targeted metagenomics approach and Illumina sequencing. As anticipated, higher variation and lower sensitivity were observed at lower coverages and allelic frequencies. We found that detection of all low-frequency variants at an abundance of 10, 5, 3, and 1%, requires at least a sequencing coverage of 250, 500, 1500, and 10,000×, respectively. Although increasing variability of estimated allelic frequencies at decreasing coverages and lower allelic frequencies was observed, its impact on reliable quantification was limited. This study provides a highly sensitive low-frequency variant detection approach, which is publicly available at https://galaxy.sciensano.be, and specific recommendations for minimum sequencing coverages to detect clade-defining mutations at certain allelic frequencies. This approach will be useful to detect and quantify low-frequency variants in both diagnostic (e.g., co-infections and quasispecies) and wastewater [e.g., multiple variants of concern (VOCs)] samples.